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Abstract 

Quantum dynamics of a charged particle in a 2D lattice subject to magnetic and electric fields is a 
rather complicated interplay between cyclotron oscillations (the case of vanishing electric field) and 
Bloch oscillations (zero magnetic field), details of which has not yet been completely understood. 
In the present work we suggest to study this problem by using cold atoms in optical lattices. We 
introduce a ID model which can be easily realized in laboratory experiments with quasi ID optical 
lattices and show that this model captures many features of the cyclotron-Bloch dynamics of the 
quantum particle in 2D square lattices. 
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I. INTRODUCTION 



A charged particle in 2D periodic potentials, which is subjected to in-plane electric field 
and normal to the plane magnetic field is a problem of lasting fundamental interest because 
of its relation to the quantum Hall effect. This problem was considered in several different 
approximations in the solid-state physics literature with the main emphasis on the energy 
spectrum of the system or, more precisely, on the electron density of states, which is a 
measurable quantity. Among these approximations the most popular is the tight-binding 
approximation, which amounts to a truncation of the Hilbert space of the single-particle 
Hamiltonian to the lowest Bloch band. In the case of zero electric field this approximation 
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results in the celebrated Hofstadter's butterfiy spectrum [Ij, which is parametrized by the 
^J'eierls phase a - number of the magnetic fiux quanta through the unit cell of the lattice 
21 . In the opposite case of zero magnetic field the spectrum is the so-called Wannier-Stark 
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ladder [3|, , which is parametrized by the angle d between the electric field vector and the 
crystallographic axis of the lattice. The case when both fields are present is more subtle 
and, to the best of our knowledge, was analyzed for the first time only in 1995 jsl. 

Complementary to the spectral problem is the wave-packet dynamics of the particle. The 
main question one addresses here is whether a localized packet remains localized in the 
course of time or it spreads over the lattice. If for electrons in a solid crystal this problem is 
of pure academic interest, it appears to be of experimental relevence for cold atoms in optical 
lattices because in this system the wave-packet dynamics can be easily tracked by taking a 
picture of the atomic cloud after a given evolution time. A recent example is Ref. iGl where 
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;he authors realized the Aubry- Andre model [7[ (which coincides with Harper's Hamiltonian 
8| for the model parameter A = 2) by loading cold atoms into the quasiperiodic ID optical 
lattice. It was confirmed that the atomic cloud remains localized for large A and spreads 
over the lattice for small A. 

Present research in cold atoms physics is also focused on the problem of generating 
synthetic magnetic fields, which could impart a ^orenz-like force to otherwise neutral atoms 
in motion (see, for example, the recent paper [9;] and references therein). Since the electric 
field for cold atoms is easily mimicked by, for example, the gravitational force, this will 
open a perspective for studying the 2D wave-packet dynamics in the Hall configuration. 



Theoretically, this problem was analyzed in much details in our recent publications lOl. 
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The main message of the present work is that many (although not all) theoretical predictions 



of Refs. 
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ll| can be verified by using the driven ID lattice instead of the 2D lattice with 
an artificial magnetic field. The proposed experiment is a modification of the laboratory 
setup [6i], where one substitutes one of the stationary lattices by a moving lattice. 

II. WAVE-PACKET DYNAMICS IN THE HALL CONFIGURATION 



To make the paper self-consistent we summarize the results of Refs. |ld . Using 
the tight-binding approximation the 2D wave-packet dynamics in the Hall configuration is 
governed by the following Schrodinger equation: 

ihipi,m = - Y (V^/+i,m + h.c.) - ^ (e*^™V«,m+i + h.c.) + ea{FJ + Fym)i)i^rn ■ (1) 

In this equation ipi^rn are the wave function probability amplitudes for the lattice site {l,m), 
Jx and Jy the hopping matrix elements along the x and y axis (in what follows we assume 

= Jy = J for simplicity), e the charge, a the lattice period, and Fy components of 
the electric field vector F, a = eBo? /he the Peierls phase, and we use the Landau gauge 
A = -6(0, X, 0) for the vector potential. 

The main conclusion of Ref. pii] is that the system ([T]) has two qualitatively different dy- 
namical regimes, depending on the inequality relation between the electric field magnitude 
and the quantity Fcr = l-naJ jea. In the strong field regime, F > Fcr, the time evolution of a 
localized wave-packet is defined by the Bloch dynamics. For vanishing magnetic field these 
would be Bloch oscillations, where the packet oscillates near its initial position with the 
Bloch frequencies u^^y = eaFr^^y/h and amplitudes proportional to J^.y/eaF^ y. The obvious 
exception from this oscillatory behavior is the case where the vector F points the x or y di- 
rection. Here the packet spreads ballistically in the orthogonal to the field direction, with the 
rate defined by the hopping matrix element. A finite magnetic field (nonzero a) 'generalizes' 
this exclusion to the cases where the vector F points rational directions, i.e., F^/Fy = r/q 



with r, q being co-prime numbers [Uj . However, now the rate of ballistic spreading in the or- 
thogonal direction is suppressed by the factor proportional to {J / eaFY^'^~^ . In practice this 
functional dependence of the suppression coefficient implies that the wave packet spreading 
can be detected only for rational directions with small prime numbers r and q. 

In the opposite limit of weak electric fields, F < Fcr, the time evolution of a localized 
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wave-packet is defined by the cyclotron dynamics. Namely, the packet moves in the orthog- 
onal to F direction with the drift velocity f *, 

v"- = ea^F/ha = Fc/B , (2) 

in close analogy with a charged particle in free space under the effect of the crossing electric 
and magnetic fields. However, the presence of the lattice imposes a restriction: this regime 
occurs only for the subspace of initial conditions spanned by the family of transporting states. 
For generic initial conditions the packet typically splits into several sub-packets moving in 
the orthogonal direction with different (both positive and negative) velocities. 



III. THE ID APPROXIMATION 



In this section we introduce the ID approximation to the above 2D problem. First we use 
the unitary transformation, '4'i,m{t) exp[—i{uJxl + ujym)t]ipi^rn{t) ■, after which the electric 
field appears as periodic driving of the system with Bloch frequencies and Uy-. 

m,m = -y (e-^'^-Vz+i.rn + h.c.) - ^ (c^^^™'"-^* + h.c.) . (3) 

Let us now assume a situation where the wave function is uniform along the y axis, i.e., 
i^i,m{t) = L~^/%i{t). This reduces ([3]) to the following ID Schrodinger equation for the 
complex amplitude hi. 

ihbi = -y (e"''^"*6i+i + h.c.) - Jy cos(27ra/ - ujyt)bi . (4) 



Although the above assumption is rather specific, it was shown in Ref. [10| that well 
approximates the integrated probability Pi{t) = J2m \4'i,m(t)\'^ also in the case of a locahzed 
2D wave packet, if its size exceeds the magnetic period d = a/a. Thus we can use the results 



of Refs. 
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ll| to explain dynamical regimes of the system (jl]) and, vice versa, to verify 
theoretical predictions of the cited papers by studying the wave-packet dynamics of this 
one-dimensional system. Note that in terms of Eq. @ the weak and strong field regimes 



correspond to slow driving, w < 



cri 



00 = ^ + uj'^ , oocr = ^iraJ/h , (5) 
and fast driving, u > Ucr, respectively. 



The system (jl]) coincides with that studied in the laboratory experiment 6| with two 
minor modifications. First, now the secondary (in terminology of Ref. {s]) lattice moves 
with the constant velocity relative to the primary lattice. This can be done by linearly 
chirping the frequencies of two counter-propagating waves which form the secondary lattice 



12| . Second, the hopping term in Eq. (jlj contains an oscillatory phase. This phase can be 
introduced by the same techniques which one uses to induce Bloch oscillations of cold atoms 
(for example, by employing the gravitational force for vertically oriented lattices). 

IV. THE REGIME OF SLOW DRIVING 

A way to check that the ID system (jlj) correctly captures the features of the 2D system 
([1]) in the weak field regime (which now assumes u < Ucr) is to propagate the transporting 
state. In the 2D lattice it moves with the drift velocity ([2]) in the orthogonal to F direction. 
For the considered system (jl]) this means that one can construct a coherent wave packet, 
which will travel across the lattice with the constant velocity 

V = aUy/2na . (6) 

Figure [T](a) shows the result of numerical simulation for the initial Gaussian wave packet , 

bi{t = 0) ~ exp(-/V2a2) , 



with the width a = {27ra^J Jy/ J^) This packet approximates the ground Wannier state 
for the potential V{1) = —Jy cos(27ra/) and is a ID analogue of the 2D transporting state 



13|. It is seen in Fig. [U^a) that the 



with the same dispersion in the orthogonal direction 
packet indeed moves with the velocity given in Eq.([6]). 

To create the narrow coherent atomic wave-packet might be a problem in a laboratory 
experiment. For this reason from now on we focus on the case of thermal atomic cloud, 
which corresponds to a wide incoherent Gaussian wave packet. We simulate the dynamics 
of this incoherent packet by assigning random phases to probability amplitudes of the initial 
Gaussian packet and averaging the result over different realizations of these random phases. 
The right panel in Fig. [1] shows the typical evolution of a wide incoherent packet. Unlike 
the case of a narrow coherent packet, now the first moment Mi = ^l^zP remains constant 



while the dispersion a{t) = J M2 — Mf grows linearly in the limit of large times. 



V. THE REGIME OF FAST DRIVING 



The characteristic feature of the cyclotron-Bloch dynamics for F > Fcr is the strong 
dependence of the rate of spreading on the field direction or, what is the same, on the ratio 
between the two Bloch frequencies. The ID system (jl]) fairly reproduces this dependence. 
As an example, in Fig. [2] we depict the coefficient A for the asymptotic linear grows of 
the wave-packet dispersion, a{t) ^ At, as a function of the frequency u for Ux/ojy = 
and LOx/oJy = 1. In the former case, the rate of ballistic spreading is seen to approach the 
constant value A = J.j./ \/2h, while in the latter case it decreases as 



where v = 1. In the general case of arbitrary rational ratio uj^/uy = rjq the exponent is 
= r + g — 1 and one has to evolve the system for algebraically large times to reach the 
asymptotic regime. Finally, for irrational oJx/ojy there is no asymptotic linear growth in a{t) 
but oscillations in time (see Fig. [3]). 

VI. FINITE EVOLUTION TIMES 

In the preceding sections we discussed asymptotic dynamics of the system. Clearly, in 
a laboratory experiment the system evolution is restricted to some maximal time interval, 
which may be not large enough to speak about the asymptotic regime. Nevertheless, all 
effects mentioned above are well observed also for short evolution times. To support this 
statement Fig. H] shows the wave-packet dispersion at t = SOTj as the function of u for three 
different ratios between driving frequencies. Two dynamical regimes, which are separated 
by the critical frequency ([5]), can be easily identified and for u > Ucr one clearly sees the 
difference between rational and irrational Ux/ujy. In addition to Fig. HJ the panels (b-d) in 
Fig. [5] depict populations of the lattice sites at the end of numerical simulations for u = 1. 
(Tiny wiggling of the curves is an artifact due to the Monte-Carlo method of simulating 
the dynamics of an incoherent packet.) We note that particular shapes of the packets seen 
in the figure are rather sensitive to variations of the system parameters and the evolution 
time. This sensitivity, however, disappears for integrated characteristics like the wave-packet 
dispersion. 




(7) 
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VII. CONCLUSION 



We have studied numerically the dynamics of non-interacting cold atoms in the driven ID 
optical lattice with a particular driving. This driving assumes the presence of a static force, 
which we characterize by the parameter - the Bloch Frequency, and a shallow secondary 
lattice with a larger period d = a/ a, which moves at constant velocity auy relative to the 
deep primary lattice. It is shown that this system well reproduces many features of the 
cyclotron-Bloch dynamics of the quantum particle in a 2D lattice. In particular, we find 
two qualitatively different dynamical regimes of the ID system, depending on the driving 



frequencies. In the first regime (slow driving, u = + ujy < Uar) a cloud of non-condensed 
atoms spreads ballistically across the lattice with a rate proportional to the frequency Uy. 
In the second regime (fast driving, u > Ucr) the size of the atomic cloud oscillates in 
time if Ux/ujy is an irrational number, while for rational u^/ujy = r/q these oscillations are 
accompanied by slow ballistic spreading with a rate inversely proportional to the frequency 
u in the power u = r + q — 1. 

We conclude the paper by a remark concerning the commensurability condition between 
the lattice periods, i.e., the rationality condition on the parameter a. This condition is 
known to play a crucial role in the case of stationary lattices Is, 7|. Unlike this situation, 
for driven lattices we did not find the commensurability condition on the lattice periods to 
be of any importance, although we do not exclude a possibility that it might be important 
in some particular situations. 
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FIG. 1: Space-time plot of the dynamics of the narrow coherent (left panel) and a wide incoherent 
(right panel) Gaussian wave-packets. Parameter are Jx = Jy = 1, a = 1/10, u = 0.1, and 
LOx/oOy = 0. The time is measured in units of the tunneling period Tj = h/ J. 
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FIG. 2: The proportionality coefficient A = A{uj) for the asymptotic hnear growth of the wave- 
packet dispersion, a{t) ^ At/Tj, for co ^/lo ,, = and uJx/^y = 1- The dashed hnes are analytical 



estimates for large u according to Ref. 
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FIG. 3: The wave packet dispersion as the function of time for ujx/^y = 1 (dash-dotted line), 
ijjxl^y = 18/19 (solid line), and ojx/^y = (\/5 — l)/4 (dashed line). The other parameters are the 
same as in Fig. 1(b) yet uj = 1. 



9 





-200 



200 




-200 



200 



FIG. 5: Population of the lattice sites at the end of numerical simulation for w = 1. Panels (b-d) 
correspond to the three cases considered in Fig. 4, the panel (a) shows the initial packet. 
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